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Abstract 

We examine phases of the Shraiman-Siggia model of lightly-doped, square 
lattice quantum antiferromagnets in a self-consistent, two- loop, interacting 
magnon analysis. We find magnetically-ordered and quantum-disordered 
phases both with and without incommensurate spin correlations. The quan- 
tum disordered phases have a pseudo-gap in the spin excitation spectrum. 
The quantum transition between the magnetically ordered and commensurate 
quantum-disordered phases is argued to have the dynamic critical exponent 
z = 1 and the same leading critical behavior as the disordering transition in 
the pure 0(3) sigma model. The relationship to experiments on the doped 
cuprates is discussed. 
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I. INTRODUCTION 



The last few years have seen numerous experiments examining the magnetic properties of 
the doped copper-oxide compounds in some detail. However, our theoretical understanding 
has not kept in step, partly due to the numerous competing effects and ensuing complexity 
of these materials: at low temperatures the effects of disorder are paramount, and at larger 
doping there is the onset of superconductivity. 

On the theoretical side, a popular model for investigating the interplay of doping and 
antiferromagnetic spin correlations has been the t-J model. Numerical and high temper- 
ature series studies of this model have been especially valuable as a testing ground for 
various theoretical ideas. In an important advance, Shraiman and SiggiaEI (SS) proposed 
a phenomenological description of the long-wavelength interplay between spin and charge 
transport in this model. There are numerous theoretical reasons for believing that their 
long-wavelength model is a correct description of the t-J model at low temperatures. It is 
expected that the SS model is quite robust, especially at low doping concentrations, and will 
describe a whole class of doped antiferromagnets, not just the t-J model. However, at suffi- 
ciently large doping, an approach based upon long-wavelength distortions of the background 
antiferromagnetic order must eventually become invalid; the mechanism of this break-down 
of the SS model is not understood and remains an important open problem. This paper 
shall examine the SS model in a new approach which is designed to explicate the nature of 
the quantum-disordered phases i.e. the phases with no long-range magnetic order. These 
phases appear at a reasonably small doping concentration, where it is expected, though not 
certain, that the mapping between the SS and t-J models is still valid. We will also discuss 
the relationship of our results to other theoretical work and some experiments in Section 

We begin by writing down a simplified version of the action S of the Shraiman-Siggia 
model: 

S = dr J d^r[Sn + Sf + S,], (1.1) 

with r the Matsubara time, r = {x,y) the spatial co-ordinates, P = l/{kBT), and T the 
absolute temperature. The first part, S'„ describes fiuctuations of the antiferromagnetic order 
parameter n^. Here is a 3-component vector and is taken to have unit length 

En,^ = l- (1-2) 

e=i 

We then have 

S,.-|((4.n.f + (W + ^(|^)^) (1.3) 

with the bare spin stiffness and Cq the bare spin-wave velocity. The momentum of the 
field is restricted to be smaller than an ultraviolet cutoff A, which is also the scale at which 
the coupling constant are defined. 

The action of the fermionic dopant holes is given by Sf. The holes are describes by 
fermionic spinor fields where a =1, | is the spin index, and v is 'valley' index. The valleys 



2 



are regions in the Brillouin zone around the minima of the fermion dispersion spectrum. For 
the square-lattice t-J model there is a great deal of evidence that there are two vallleys 
centered at the points (vr/2, ±vr/2) in the reduced Brillouin zonei'i. We will rotate our co- 
ordinate system from the conventional one, so that our x-axis is at an angle of 45 degrees 
to the axes of the square lattice, and the principal axes of the valleys are along the new x 
and y axes; the index v will therefore take the values v = x,y. Further, we will measure 
momenta from the center of the valleys. With these conventions, the fermionic action Sf is 

5f = $L ( — - — - -TT^ 

2mh oy^ / 




+ $t (14) 
^ 2mhdx^ 2midy^) ^ ' 

where mi, nih are the light and heavy masses of the hole. 

It is important to realize that there is no simple relationship between the bare fermionic 
field cJq, of the t-J model {i is the site index) and the continuum fields (n^, $^,a) of the 
hydrodynamic SS model. Crudely speaking, one may consider this field transmutation as 
a form of spin-charge separation in which the spin-1/2, charge e fermionic field has 
separated into the spin-1, charge 0, bosonic rie quanta and the spin-1/2, charge e, 
fermions. 

Finally, we will consider only the leading term in the coupling, Sc between the holes and 
the Hi field, which is responsible for inducing local incommensurate spin correlations. It is 
expected that the remaining terms in the SS model are innocuous and do not change the 
results of this paper qualitatively. We have 

Sc = K ($L^Q/3'^^/3) {(^empUmd^np) (1.5) 

where are three Pauli matrices, and k is a coupling constant. A crucial feature of the 
SS model is that there is no three-body coupling like rii^l^a^p^^is between the fermions 
and n£ quanta!: such a term is forbidden by a sublattice interchange symmetry of the 
SS mo delS under which n£ change sign while the remain invariant. Other models of 
nearly antiferromagnetic fermi liquidsji posses such a term and as a result, have quantum 
disordered phases with rather different properties. 

We wish to clarify a crucial point about our particular form of S at the outset. SS have 
shown that there is a mapping, in principle exact, to an alternative form of S in which the 
degrees of freedom are a spin-1/2, charge 0, complex scalar Za with a =t, i and spinless, 
charge e fermions "^^a where a = A, B is a sublattice index. The (n^, $^a) fields are related 
to the {za, "^va) fields by 

"'->v- (:::)^ (;■;/)(*::) 

Note that there is a U{1) gauge transformation on the (zaj'^va) fields which leaves the 
(ui, ^va) fields invariant: 



■ ^f.A^'^vAe-'^ ; ^.Bc't (1.7) 



3 



where has an arbitrary dependence on spacetime. The action in terms of the {za, "^va) 
will therefore also have to be invariant under the gauge transformation; there is no such 
restriction on the action in the (n^, ^^a) variables. This latter absence of gauge restrictions, 
and associated long-range gauge forces will be quite useful to us in our analysis. 

SS go on to state that equivalent results are obtained in computations using either 
choice of fields. In principle, this statement is correct. In practice, however, one is usually 
restricted in the analysis to perturbation theory, in which the quantum numbers of the low- 
lying excitations are essentially identical to those of the degrees of freedom in the action. 
Equivalent results in both theories have been obtained in the vicinity of magnetically ordered 
states, where one is performing a small fluctuation, spin-wave analysis. On the other hand, 
the choice of fields is expected to have dramatic consequences in a quantum-disordered state. 
Kane et. al. and other^ made a choice of fields equivalent to the (zq, formulation of 
the SS model. Thus, not surprisingly, their Schwinger-boson mean-field theory yielded a 
quantum disordered phase with deconfined, spin-1/2, charge 0, bosonic spinous (the Za) and 
spinless, charge e fermions {^va)- 

In this paper, we use another approach and shall examine the phases of the SS model 
which are obtained naturally in the ira^, ^va) formulation. Some of the results of our cal- 
culations were noted some time agoB. While, in the end, we have no formal justification 
for claiming that the (ra^, ^^a) formulation is more accurate than the {za-, '^va) approach, 
we can offer the following motivations. The field formulation has been quite successful 
in describing the undoped, frustrated antiferromagnet. The long wavelength action of the 
undoped antiferromagnet! is simply the 0(3) sigma model, Sn-, and it is expected to display 
a quantum-disordered phase in which massive n quanta form the lowest excitations and 
carry spin 10. Consistent with this, recent investigations of the quantum-critical behavior 
in these systems have argued for the superiority of the n-field based approach, and have 
successfully explained a number of experimental and numerical computations on the square 
lattice antiferromagneti0. Further, careful analysis of fiuctuationJli in the Zo-based the- 
ories of the undoped antiferromagnet has yielded quantum disordered phases in which the 
quantum numbers of the low-lying excitations are identical to those n-field, which means 
that fluctuations in fact confine bosonic spinous into 5* = 1 particles. It is then natural 
to explore the consequences of doping in a model in which the correct physics in the limit 
of zero doping is captured most directly i.e. in the (n^, approach. We shall argue 
later in this paper, that the results of such an investigation are consistent with the available 
numerical and experimental data on doped antiferromagnets. 



A. Summary of Results 

We now discuss the main results of our calculations. We will distinguish the quantum 
phases by the properties of the equal-time ni correlator 

~S{q) = {W{q)?) (1-8) 

as T — >■ 0, as a function of the momentum q. We emphasize that S{q) is not the full structure 
factor S{q) measured in neutron scattering experiments; S{q) will contain additional terms 
involving the contribution of the spin-1/2 <^ya fields. Phases with magnetic long-range 
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order have a delta function term in S{q) at T = 0; this delta function is at g = in the 
commensurate long-range-ordered phase (hereafter referred to as CLRO) and at g 7^ for the 
case of incommensurate long-range-order (ILRO) (g is measured from (vr, vr)). The quantum 
disordered phases have no delta-function terms at T = 0, but only a peak of finite, though 
possibly small, width. This peak is at g = in the commensurate quantum-disordered phase 
(CQD) and at g 7^ in the incommensurate case (IQD). In agreement with the SS analysi^l 
and experimental results^, the peak in the ILRO and IQD phases was found to occur along 
the conventional (1, 0) and (0, 1) axes of the square lattice (these are the (1, 1) and (1, —1) 
axes in our rotated co-ordinate system). 

We will have little to add here to existing studiellEHiel 

of the properties of the 

magnetically-ordered phases (CLRO and IRLO): their low-lying excitations are spin- waves 
involving long-wavelength deformations of the ordered state. Our focus will mainly be upon 
the new quantum-disordered phases (CQD and IQD) and their unusual properties. The 
quanta in both phases were found to be fully gapped. The low-lying excitations in the 
sector consist of a triply-degenerate spin-1 particle with a finite energy. However, the spin- 
1/2 particles continue to form a Fermi sea which possesses gapless fermionic excitations 
with charge e and spin-1/2. Despite the presence of these gapless excitations, the gap 
is robust as there is no term in the SS model which permits the decay of a ni quantum to 
a fermion particle-hole pair. The importance of the absence of the three-body term noted 
above, is now evident. Taken as a whole, the model thus only has a pseudo-gap to spin 
excitations in the CQD and IQD phases. One of the consequences of the presence of the 
gapless spin-1/2 fermions is that the uniform spin susceptibility of the CQD and IQD phases 
will be finite at T = due to the Pauli contributionQ. We also note, that our calculation has 
completely neglected the effect of Berry phases; in the context of undoped antiferromagnets 
it has been arg ued0 that Berry phases should induce spin-Peierls long-range order in the 
CQD phases. It is possible that such spin-Peierls order will also exist in the CQD phase of 
the doped antiferromagnet. 

We undertook a partial numerical survey of the phase diagram of the SS model as a 
function of k and the hole density. Parameters were always chosen so that the zero 
doping state was CLRO. This CLRO state was always found to be stable over a small, but 
finite, doping concentration. Over some of the regime examined, the sequence of phases 
with increasing doping was CLRO - CQD - IQD. We studied the T = quantum transition 
between the CLRO and CQD phases and will present evidence indicating that it has dynamic 
critical exponent z = 1 and the same leading critical behavior as the transition in the pure 
0(3) sigma model; however, the corrections to scaling in the two models were found to 
be quite different. The boundary between the CQD and IQD phases is an example of 
a disorder lineS: our calculation only found a non-analyticity in the dependence of the 
structure factor on the bare coupling constants at the disorder line, but no strong long- 
wavelength fluctuations. 

In a region of the phase diagram with k large, we found the sequence CLRO - ILRO 
with increasing doping. In principle, there should eventually be a ILRO to IQD transition, 
but for the parameters examined, we did not find one before a doping level where the 
incommensuration wavevector was almost as large as the momentum upper cutoff. 

We will argue from the above numerical results, and from theoretical considerations, that 
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there is a Lifshitz point in the k plane where all the four phases - CLRO, ILRO, CQD, 
IQD - meet. Some properties of this multicritical point will be discussed. 

We also have extensive results on the temperature dependence of equal-time correlation 
functions in the various phases. In particular, the temperature dependences in the spin 
correlation length and the structure factor are quite instructive, and will be described later. 



II. CALCULATIONS 

For the case of the undoped antiferromagnet, the expansion on the 0{N) non- 
linear sigma model offers a convenient and accurate method for exploring properties in the 
vicinity of T = quantum transition JHl. The extension of the 1/N expansion to the doped 
antiferromagnet is however not straightforward because of the presence of the third-rank e£mp 
tensor in Sc (Eqn. (|1.5|) ), which is special to the case = 3. Even with this complication, it 
is still possible to justify perturbative 1/N calculations, although in a rather inelegant way: 
after the fermions have been integrated out, the coupling k in the effective action of the n 
field has be to scaled by l/N^^ where < /i < 1/2. Not much is learned from this extension 
to general A^, and we will therefore spare the reader the details. We will be satisfied, 
instead, in restricting our discussion to the special case of = 3, and viewing our 1/A^ 
calculation as a physically motivated, self-consistent, interacting magnon approximation. 
The magnon-magnon interactions are computed in a manner which is directly inspired by 
the 1/N expansion of the undoped antiferromag net§. 

An important property of our approach is that spin-rotation invariance is explicitly 
preserved at all stages of the calculation. This is crucial for a proper study of the quantum 
disordered phases of the model, especially when the n quanta acquire a gap. We thus 
expect our approximations to work best in the quantum-disordered phase, at the T = 
quantum transition, and in the intermediate temperature quantum-critical regioni'Q. At 
the same time, the low temperature properties in a region with magnetic long-range order 
(the renormalized classical regionll) in the ground state, may not be well described. Even 
in the undoped antiferromagnet, the 1/A^ expansion is singular in the renormalized-classical 
region, and a careful interpretation of the results is requiredliil. We therefore will focus below 
here will be mainly on the quantum-disordered phases phases. 

We begin by expressing the action in suitable dimensionless parameters. We rescale 
lengths such that the upper-cutoff in momentum space for the n field is 1; thus 

r 

r-^- (2.1) 
Similarly, the times r are rescaled so that the bare spin-wave velocity of the n field is unity: 



After rescaling the temperature 



T^T^ (2.3) 
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we then have the following modified form of 5'„ 

Sn = ^ Ud^ne)' + (dyne)' + (S.n,)^) (2.4) 
where the dimensionless coupling constant g is given by 

« - ^ (2.5) 

The fermionic action will retain its form after rescaling the field $ A$ and rescaling to 
the dimensionless effective masses 

mh,i mh,i — (2.6) 
Co 

Finally in Sc we replace k — > kcq. 

We will impose the constraint ( |1.2| ) by a Lagrange multiplier field A. Thus we need to 
evaluate the functional integral 



Z = J VniV<^ya'DXexp(^-S 



-I 



dr / d''x—{nl - 1) (2.7) 



J 2g 



It is now possible to set up a rotationally-invariant, diagrammatic expansion of all ob- 
servables associated with Z. We will work at finite T, and so no breaking of spin rotation 
invariance can occur; the properties of the ground state will be elucidated by taking the 
T — > limit. We treat the A field in much the same way as in the undoped systen£3. We 
assume that fluctuations of iX occur about a saddle point value A; we therefore write 

iX = X + iX (2.8) 

where A is the fluctuating part of A. The value of A is to be determined at the end of the 
calculation to satisfy the constraint ( |1.2| ). The diagrammatic expansion now has three bare 
propagators: the conventional Green's function of the fermions ^ya, the propagator, of 
the rii field 

G%q,tu^)= . (2.9) 



{q is the wavevector and u;„ is a Matsubara frequency) and the propagator l/II of A, witht2l 

U{q,iu;n)=Tj2 J ^^"(A; + g, 2e„ + za;„)G°(A;, ^e„). (2.10) 

There are two interaction vertices: the four-body (n — n — $ — $) coupling in Sc and a three- 
body (X—n—n) vertex with the value i/ {2g). Finally, there is a rule to prevent over-counting: 
no A propagator can be followed by a bubble consisting just of two G° propagators. 
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We may now write the fully-renormalized correlator of the n field in the form 

G{q,iuJn) = — — ^—^7 — : — ^ , . (2.11) 

where E is the self energy and the 'mass' m is given by 

m2 = A + S(0,0) (2.12) 

The lowest order contributions to S from magnon-magnon (S„) and magnon-fermion (Sj) 
interactions are shown in Fig ^ and their values are 

S = S„ + 

(fk G'^{k + g, ien + iujn) 



47r2 n(A;,ze„ 



En{q,iuJn) = ^Tj2j 

— (2g^ + k^fxv{k,ien)G^{k + g,2e„ + iun) (2.13) 

where Xv is polarization of the fermion in valley v. We used the following expression, 
appropriate for an elliptical Fermi surface: 



Xx{q,ien) = Xo{T) 



1 N 

\ 1 1 /2 



(2.14) 



and similarly for Xy The Fermi velocity Vp-, the Fermi wavevector kp and the doping 
concentration 5 are related by the equations 

vp = ^^ ■ 5 = ^. (2.15) 

and the polarization is taken to vanish unless 

(g^(mjm0^/' + gJ(^,/m^)i/2)i/=^ < 2fc^ (2.16) 
The prefactor xo we chose as the compressibility of a free Fermi gas at a temperature T: 

^.(T) = (l - e-'=-^-/(2^)) (2.17) 

Our approach consisted of using the above approximation for S and then solving the 
constraint equation (|1.2| ), or 

3^?^E / ^G'(g,z^„) = l (2.18) 

for the value of (or equivalently A). The dependence of S on G"^ was made partially 
self consistent by replacing A by in (|2.9|), thus using G'^{q,iun) = 1/(9^ + "-^n + 
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in (|2.1CI| ) and ( |2.13|) . A fully self-consistent approach would require we replace G'^ by G 
in these equations; this is computationally much more difficult and was not numerically 
implemented. Our approximation thus amounts to replacing G° by G, but then ignoring 
the momentum and frequency dependence of the self energy in G. For the most part, this 
omission is not expected to be serious, as corrections can be organized order by order in k. 
However, we cannot rule out the possibility, especially in the quantum-disordered phases, 
that there is some entirely different, possibly gapless, solution of the fully self-consistent 
equations; such a solution will clearly be non-perturbative in k. We also note here that in 
our analytical considerations below of the boundaries between the phases, we will include 
the full G in (57X31). 



The numerical determination as a function of T, S was carried out on a HP- RISC 
workstation. A meaningful solution always existed at all finite T, with no phase transitions 
as a function of T or 6. Phase transitions are however present at T = 0, and were examined 
by studying the T — > limit of our solutions. The computations required about 3 weeks of 
computer time. 



III. RESULTS 

We now describe the results of our numerical calculations. The nature of the ground 
state can be determined from the values and T dependences of qc and m where q = qcis the 
location of the maximum of G{q,iuJn = 0), and 

= + E(g„icj„ = 0) - E(0,0) (3.1) 

It is easy to see from (|2.11| ) that m is roughly the inverse correlation length ('roughly' 
because this neglects dH/dq'^; including this term yields corrections of the order of unity 
which are not strongly T dependent). The various phases can be identified by studying the T 
dependence of m as T — >■ 0, as will be described below. The values of qc were approximately 
T independent and distinguish between commensurate and incommensurate phases. 

Two samples of our results are contained in Figs. |^ and § which plot the T and doping 
dependence of m for two sets of coupling constants. For completeness we also show in Figs |^ 
and ^ the values of m for the same samples. We will now describe the properties of the 
phases in these figures and follow that up with some general discussion on the nature of the 
quantum transitions between them. 



A. Long-range-ordered states 

The states with magnetic order are expected to have m — as T 0. In particular, 
the low T dependence should beEl 

M^exp(-^) (3.2) 

where ps is the fully renormalized spin stiffness. Numerical solutions at very lower T took 
longer times to converge, so it was difficult to see this exponential behavior in some of the 
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doped samples. We simply identified the samples in which m vanished with an upward 
curvature as T — 0, as possessing magnetic long-range order. Further states with = 
(g, ^ 0) were identified as CLRO (ILRO). 



B. Quantum disordered states 

These states have fn saturating at a finite value as T — > 0, which is roughly the gap, 
A ~ m(T ^0), in the n sector. Again, the value of Qc distinguishes between the CQD and 
IQD states. 

We examine the nature of the spin correlations at a point in the IQD phase by plotting 
the n field contribution to the structure factor, S{q), 

~S{q)=gTY,G{q,u^) (3.3) 

in Fig. p. Notice that there is strong overlap between the peaks at high temperature. 
Upon lowering the temperature, the peaks first sharpen considerably, but then their width 
saturates. 

Let us discuss the form of the n spectrum at T = in the CQD phase. We will focus on 
real frequencies, uj just above the gap A, and small momenta q. The magnon contribution 
to the self energy, in (|2.13|) does not acquire an imaginary part until uj = 3A and can 



therefore be completely ignored. The damping from the fermion particle-hole continuum, 
S J is however not so innocuous. We findEl 



ImEj(g, u) 



< < A 

\q\{uJ-Af 0<^-A«^ (34) 

Ai/2(^_A)5/2 ^_A»4- 
^ ^ 2A 



and Imlly(g, —u) = — ImSj(g, tu). The n spectral weight is then given by 

Giq,uj) - .2 I 2 \ , t V (3-5) 

From the above results it follows that at g = we have 

ImG(g = 0, u; > A) = ^^(c. - A) + ^(c. - A) (3.6) 

for some constants ai,a2- Thus there is a sharp spin-1 quasiparticle peak, and a second 
background term which is a direct consequence of the coupling of the n quanta to the 
particle-hole continuum. At small, but finite g, the sharp peak moves too; ~ A + gV(2A) 
and acquires a finite width; there is absorption at all frequencies greater than A. 

The spectral properties of the IQD phase are essentially identical except that the role of 
the point g = is replaced by g = q^, in obtaining this result it is, of course, necessary to 
replace G° by G in (pA3|) . 
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C. Quantum transitions 



Our results in Figs. |^ and ^ show two sequences of quantum transitions with increasing 
doping: CLRO-CQD-IQD and CLRO-ILRO. In the second case there should eventually be 
a ILRO to IQD transition, but for the parameters examined, we did not find one before a 
doping level where the incommensuration wavevector was almost as large as the momentum 
upper cutoff. We now present a theoretical analysis of some issues raised by the existence 
of these quantum transitions. 



1. CLRO to CQD quantum transition 

An important ingredient in determining the universality class of this transition is the 
analytic structure of the n field self-energy S j as T ^ in the Neel phase and the quantum- 
critical point. We consider ( p.l3| ) in the limit m — 0, and T — when the frequencies 
become continuous variables, and the Matsubara summations can be converted to integra- 
tions. It is evident that T,f{q,iuj) is an even function of uj. Moreover, it is not difficult to 
show that there are no infra-red divergences in either 5Sj/(9g^|g=^=o or dT, j / duj'^\g=o^^\o. 
This implies that for q, uj small we have 

^f{q,iuj) = Sj(0,0) + 6ig' + 62C^2 + ... (3.7) 

Thus the gapless fermion particle-hole sea has not induced any non-analyticities in to this 
order. There are indeed non-analytic terms present at higher order in which are signaled 
by infra-red divergences in higher derivatives of Ej; we will discuss the form of such terms 
below. For our purposes, it is sufficient to note here that all such higher gradient terms 
are expected to be irrelevant at the CLRO to CQD transition. Thus the gapless fermion 
particle- hole excitations have had a relatively innocuous effect: they have mainly lead to 
renormalizations of the spin-wave velocity and spin stiffness. The universality class of the 
CLRO-CQD transition is thus expected to be the same as that in the undoped sigma model. 
This is a transition with dynamic critical exponent z = 1 and its leading universal properties 
have been discussed in some detail by Chubukov et. a/.0. All of the scaling functions of 
Chubukov et. a/.0 should therefore also apply to the present doped antiferromagnet. The 
main effect of the fermions has been to change the value of the effective coupling constant 
and renormalize the spin-wave velocity. Consistent with this identification, observe the linear 
dependence of m with T in Fig. ^ at fci? = 0.2 over a wide temperature region. This value of 
kp places the system quite close to the quantum-critical point as the value of A is very small. 
At the quantum-critical point of the sigma model, it is predicted that = CqksT /he with 
Cq f» 1 a universal number. The slope of m versus T &i kp = 0.2 in Fig. ^ is about 0.65 
- this matches with the expected result if there is renormalization of spin-wave velocity 
c/cq ~ 1/0.65. A renormalization of the spin-wave velocity of order unity is to be expected, 
as the fermionic polarization x in ( |2.17D is not suppressed by a factor that vanishes as 5 — > 0. 



Differences between the quantum transition in the doped and undoped antiferromagnet 
do however show up at the correction to scaling level. The higher-order non-analytic terms 
in Ej- will have a form which is quite specific to the doped model. One such term can be 
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obtained by analytically continuing to real frequencies and computing ImE/ at the critical 
point. We find 

ImSj(g, u;) ~ gu;^ at the quantum-critical point (3.8) 



2. Lifshitz point 

The existence of a direct CLRO to ILRO transition in Fig. ^ has a strong consequence for 
the phase diagram of the SS model. As both phases can be transformed into their quantum- 
disordered partners simply by increasing the value of g, we conclude that there must be a 
point in the phase diagram where all the four phases - CLRO, CQD, IRLO, and IQD - meet. 
Such a point is called a Lifshitz pointS. Lifshitz points have so far been studied primarily 
in the context of thermal transitions in classical spin systemJ^ll. An important result is 
that such points can exist only above a lower critical dimension determined as foUow^Hi: a 
system in D dimensions, with incommensurate instabilities in m of those dimensions, has 
lower critical dimension 2 + m/2. This result appears to be in conflict with our results here 
for the doped antiferromagnet. For we have incommensuration in m = 2 spatial dimensions, 
giving a lower critical dimensionality equal to the spacetime dimension D = 3. So how can 
a Lifshitz point exist ? 

The answer to this apparent inconsistency lies in the form of S/. At the Lifshitz point 
we clearly have 9G^^/9g^|g=o,w=o = 0. Thus the leading q dependence of at small q will 
come from higher-order terms in Sj. Let us assume that ~ cj^-t-g^ at the Lifshitz point. 
Inserting this fully renormalized G in the result (|2.13|) for S/ we find by power counting 



S/(g, 0) = S/(0, 0) + hq' + hq^-^'^ + ... (3.9) 

Consistency now demands that p = 4 — ^/2 which yields p = 8/3. This differs from the 
value p = 4 use in classical spin systems". With this modified form of G we may repeat 
the calculation of Crest and SakH and verify that spacetime dimension D = 3 is above the 
lower critical dimension which is D = 7/3. Thus it is possible to have a Lifshitz point in 
D = 3. 

Finally we note that a point where CLRO, ILRO, CQD, and IQD phases meet was also 
found in the large N, Sp{N) theory of frustrated, two-dimensional quantum Heisenberg 
antiferromagnetsEa. The nature of this point appears to be quite different from the Lifshitz 
point in the present theory. In particular, the IQD phase of the Sp{N) frustrated antiferro- 
magnetS contains deconfined, bosonic, spin-1/2 spinous, while here we have found massive, 
triply-degenerate n quanta. Perhaps related to this difference is the fact that the large 
Sp{N) theoryil finds no softening in the spinon spectrum at the CQD-IQD boundary. In- 
stead, at the boundary the parabolic spinon spectrum splits into two parabola with minima 
at incommensurate points; the curvature at the minima of the parabola always remains 
finite. Contrast this with the behavior of the n spectrum found here: the curvature at the 
minimum of the n spectrum vanishes at the CQD-IQD boundary leading to double minima 
at incommensurate points in the IQD phase. 
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IV. CONCLUSIONS 



The most notable features of our results on lightly- doped antiferromagnets are the quan- 
tum disordered phases with a spin pseudo-gap. These phases posses fully-gapped, triply- 
degenerate, spin 1 magnons, and gapless, spin-1/2, charge e fermions. Several other in- 
vestigatorsB^t^ have also recently explored models of the normal state of the lightly-doped 
cuprates which have spin gaps/pseudo-gaps. This interest in pseudo-gaps is of course mo- 
tivated by numerous experiments on the underdoped cuprates showing gap-like features 
in the normal state.&Ii There are also interesting trends in the doping and temperature 
dependence of uniform spin susceptibility of the cupratesS. 

For completeness we review some of the previous theoretical results, and point out the 
differences to our results. A number of the models^ are related to resonating-valence- 
bond type mean-field theoriesil; there is a BCS-like pairing of spin-1/2, neutral spinous 
in the normal state at low-doping, leading to a gap-like feature in the spectrum at finite 
temperature. However, unlike our results, this state extrapolates to a true gap at T = 0. 
An extension of these models to a three-band, Cu02 layer modelil did posses gapless, spin- 
1/2, charge e fermionic excitations on the oxygen sites. However the above-gap spectrum 
in all of these modelJ^^S consists of unbound spin-1/2, neutral fermions; in contrast, the 
spin spectral weight above the gap of our model is dominated by a spin-1, bosonic magnon. 
Millis and Monien have attributed the gaps to interlayer couplings in the Yttrium based 
cupratesii. The recent work of Sokol, Pines and collaborators0 is perhaps closest in spirit 
to ours, although their scenario for the mixing between the n quanta and the fermions 
appears to be somewhat different. 

Also relevant to our result is the recent high temperature series analysis of the CLRO to 
CQD transition in the t — J mo deli. This work provides some evidence in support of our 
result that z = 1 at this transition. 

Finally, if our model is to provide a complete picture of the cuprates, it should also explain 
the nature of the photo-emission spectrumE^I. For this one needs to understand better the 
connection between the bare electron ejected in the photo-emission and the (n^, <^>^,a) fields. 
This problem is currently under investigation. 
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FIGURES 



FIG. 1. Feynman diagrams for the n self-energy. The thick, full line is the n propagator, the 
thin, full lines are the hole fermions, and the dashed line is the A propagator. 



FIG. 2. Values of fn (defined in Eqn. (3J.); it is roughly the inverse correlation length) as a 
function of T for various doping levels specified by kp- The coupling constants were g = 3, k = 1.5, 
{mumhY^'^ = 0.6, mi/mh = 0.1296. The maxima in the structure factor are at g = qc{H, ibl)/-\/2 
in our rotated co-ordinate system; this places them at (vr it Qc, tt) and (vr, vr it Qc) in the conven- 
tional Brillouin zone of the square lattice. All parameters are measured in the dimensionless units 
described in Section |ll[ The small T behavior of m, and the value of Qc identifies the nature of the 
ground state (CLRO, ILRO, CQD, IQD) which is also noted. 

FIG. 3. The values of m associated with the results for fn in Fig. |2[ For the commensurate 
states m = m. 



FIG. 4. As in Fig. | but for g = 6, k = 2, (m^m/i)!/^ = 0.6, mi/rrih = 0.1296. 
FIG. 5. The values of m associated with the results for m in Fig. ^. 

FIG. 6. Contribution of the n field to the structure factor, S{q), in the IQD phase. We use 
the coupling constants of Fig. Q at fc^? = 0.3. The temperatures are (a) T = 0.5, {h) T = 0.2, and 
(c) T = 0.02. The jagged double-peaks in (c) are an artifact of the plotting routine. 
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